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Low- frequency spin dynamics in the S — 1/2 antiferromagnetic spin-chain compound a- 
VO(POa)2 has been studied by means of 31 P NMR. The nuclear spin-lattice relaxation rate 
1/Ti at the P site exhibits H~ x ' 2 dependence on the applied magnetic field (H) at temperatures 
(T) well above the intrachain coupling strength J/kn = 3.50 K indicating one-dimensional dif- 
fusive spin dynamics. The diffusive contribution to 1 /T\ decreases on cooling as electronic spins 
acquire short-range antiferromagnetic correlations within the chain, and vanishes almost entirely 
around T ~ J/kn- This is accompanied by an apparent increase of the spin-diffusion constant 
from the value expected in the classical limit. On the other hand, the field-independent part of 
1/Ti increases with decreasing temperature, which may be a precursor for the true long-range 
antiferromagnetic ordering found below Tn = 1.93 ± 0.01 K. 
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§1. Introduction 

There has been a continued interest in the dynamics 
of low-dimensional quantum antiferromagnets at finite 
temperatures. One of the issues which has attracted 
renewed attention is a problem of spin diffusion in the 
one-dimensional (ID) Heisenberg spin chain. It is argued 
from phenomenology that at high enough temperatures 
and at long times, the spin autocorrelation function of 
the ID Heisenberg spin chain has a diffusive form 



{Si(t)Si(0))<xt-V*, 



(1.1) 



leading to divergence of the spectral density at low fre- 
quencies as bj^ 1 / 2 . However, because the diffusive form 
(1.1) is not derived from the microscopic Hamiltonian 
but is a consequence of the hydrodynamical assumption 
for the spin-spin correlation,Ef the question on the exis- 
tence of spin diffusion in ID spin chains has been stud- 
ied intensively from both theoretical and experimental 
viewpoints. One of the best studied theoretical models 
is tbGi-S-i=i4/2 XXZ chain represented by the Hamilto- 



H 



qy qV 



AS?S? +1 ). (1.2) 



At the isotropic point (A = 1) and at low enough tem- 
peratures (T <C J/fce), analytical expressions for the dy- 
namical susceptibility x(g, uS) have been derived for q ss 
and 7T, and are shown to have no diffusive pole at low fre- 
quencies^ At modestly high, or much higher tempera- 
tures compared with the intrachain coupling strength J, 
the problem is still controversial. Although the absence 
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of diffusive excitations seems to be settled for the XXZ 
chain with planar anisotropy (0 < A < 1), a definite 
answer for the absence (or presence) of spin diffusion at 
the isotropic point has not yet been given. 

Diffusive behavior of the spin-spin correlation has been 
observed experimentally in several ID spin chains via 
the w -1 / 2 resonance-frequency dependence of the nu- 



clear spin 
peratures. 



Stice rpla.Ya.tinn rate 1/Ti at elevated tem- 
@BB§ As to the 5=1/2 chain, an 
analysis of the u> dependence of 1/Ti at the Cu site 
in Sr2Cu03 based on the classical spin-diffusion the- 
ory gives an unusually high value of the spin diffusion 
constant compared with the classical value.EiP This sug- 
gests the absence of spin diffusion in the low-temperature 
limit in consistency with the field-theoreticalEr and 
perturbativair approaches. On the other hand, an 
evidence for diffusive spin transport is found from 
the proton NMR measurements in CuCl2-2NC5H5 and 
Cu(C 6 H 5 COO)2-3H20 at temperatures well above the 
intrachain exchange interactions.liiP Therefore, there 
seems to be a gap between the high- and low-temperature 
behavior of the spin transport in the S — 1/2 chains, 
necessitating further experiments on the dynamics, es- 
pecially on the temperature-dependent properties. 

Linear chains of a V 4+ ion in the compound a- 
VO(POa)2 may be a model system of an S = 1/2 Heisen- 
berg spin chain. a-VO(P03)2 belongs to the monoclinic 
space group C2/c and has the room-temperature lattice 
parameters; a =^5.140 A, b = 4.195 A, c = 9.573 A 
and f3 = 120.54° H> In the a-VO(P0 3 ) 2 structure (Fig. 
1), VO5 pyramids are stacked along the b axis to make 
up a linear chain of V atoms with the nearest-neighbor 
V-V distance of 4.915 A. The linear chains are well sep- 
arated by distorted PO4 tetrahedra in the a and c di- 
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oclinic b axis. The dark- and light-shaded polyhedra represent 

VO5 pyramids and PO4 tetrahedra, respectively. T 



rections, so that good one dimensionality in the b direc- 
tion is expected. In this paper, we report on the exper- 
imental study of the low-frequency spin dynamics in a- 
VO(P0 3 ) 2 via 31 P NMR. The relatively small intrachain 
coupling (J jk-Q = 3.50 K estimated from the present sus- 
ceptibility measurement) in a-VO(P03)2 enables us an 
experimental access to a wide range of temperatures not 
only T ~ J/ks where short-range antiferromagnetic cor- 
relations are important, but also T 3> J/k B where elec- 
tronic spins behave almost paramagnetically. Crossover 
of the dynamics between the two temperature regimes 
can also be elucidated by examining the temperature- 
dependent behavior of, for example, the nuclear spin- 
lattice relaxation rate. This type of the experiment can- 
not be done in the canonical S — 1/2 Hcisenberg anti- 
ferromagnetic spin-chain compound Sr2Cu03 which has 
a huge intrachain exchange (J/fce = 2200 K)^ being 
suitable for the study of low-temperature dynamics, and 
will give complemental information for thorough under- 
standing of low-energy spin excitations in the S = 1/2 
Heisenberg spin chain. 

§2. Experiments 

Polycrystalline samples of a-VO(P03)2 were prepared 
by a solid-state reaction method. Equimolar mixture of 
(VO)2P207 and P2O5 was fired in an evacuated silica 
tube at 750 °C for 2 days and at 900 °C for 2 days with 
the intermediate grinding. (VO)2P207 was prepared as 
described in ref. [IT]. The obtained samples were exam- 
ined by the X-ray diffraction measurement and were con- 
firmed to be a single phase. Magnetic susceptibility was 
measured using a SQUID magnetometer (Quantum De- 
sign MPMS-5s) at 0.1 T. NMR measurements were per- 
formed with a standard phase-coherent pulsed spectrom- 
eter. 31 P NMR spectrum was taken by recording the 
spin-echo signal with a Box-car averager at a fixed fre- 
quency while sweeping the external magnetic field. The 
nuclear spin-lattice relaxation rate of 31 P was measured 
by a saturation recovery method with a single saturation 
rf pulse. The measured nuclear-magnetization recoveries 
were single exponential as expected for nuclei with spin 
1/2, so that the nuclear spin-lattice relaxation time T\ 



Fig. 2. Temperature dependence of the magnetic susceptibility of 
o-VO(P03)2- The inner core diamagnetic susceptibility (—9.2 X 
10 — 5 emu/mol) was subtracted from the raw data. The inset is 
the expanded view at low temperatures. The solid line in the 
inset is a fit of the data to the susceptibility of an S = 1/2 
Heisenberg antiferromagnetic linear chain. 



can uniquely be determined as the time constant of the 
recovery curve. 

§3. Results and Analysis 

3.1 Magnetic Susceptibility 

Figure 2 shows the temperature (T) dependence of 
the bulk magnetic susceptibility (\) of a-VO(PC"3)2- 
While showing Curie- Weiss-like behavior at high tem- 
peratures, the susceptibility takes a rounded maximum 
around T max w 4.4 K as a sign of low-dimensional nature 
of the exchange coupling. The susceptibility can nicely 
be fitted to the Bonner-Fisher curve for the S =1/2 
Heisenberg antiferromagnetic spin chain (HAFC).E3 We 
estimated the intrachain coupling J between V 4+ spins 
by fitting the T dependence of the bulk x to the for- 
mulaJiS 1 



X =Xo 



k B T 

0.25 + 0.14995 x + 0.30094 x 2 
1 + 1.9862 x + 0.68854 x 2 + 6.0626 x 3 ' 



(3.1) 



Here xo is a T-independent part of x, Na is Avogadro's 
number, g is Lande's g factor, ji B is the Bohr magne- 
ton and x = J/ksT. From the least-squares fit of all 
the available data points (1.8 to 300 K), we obtained 
Xo = 1.05(8) x 10~ 5 emu/mol, g = 1.978(1) and J/k B = 
3.50(1) K. The result of the fit is shown in the inset of 
Fig. 2. xo m ay be interpreted as the Van-Vleck orbital 
paramagnetic susceptibility xvv, and the obtained value 
of xo is in a reasonable range for xvv of a V 4+ ion. The 
bulk x deviates slightly from the Bonner-Fisher curve 
below about 2.2 K which may be attributed to the effect 
of the interchain coupling. As expected, a transition to 
the long-range ordered state is found to occur at 1.93 K. 
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Fig. 4. Temperature dependence of the principal values of the 
Knight-shift tensor, (open circles) and K± (solid circles), at 
the P site in «-VO(P03)2- The inset shows the K-\ plot for each 
principal component. The solid lines in the inset are the linear 
fit of the data, the slopes of which yield the hyperfine coupling 
constants An and A i . 



Fig. 3. Temperature variation of the field-swept NMR spectrum 
of 31 P in a-VO(P0 3 ) 2 taken at 48.5 MHz. The dotted line 
indicates a zero-shift position for 31 P. The inset shows the tem- 
perature dependence of FWHM of the spectrum. 



3.2 NMR spectrum and the Knight shift 

Typical examples of the field-swept 31 P NMR spec- 
trum are shown in Fig. 3. Above 2.0 K the spectrum 
exhibits an asymmetric pattern resulting from an axially- 
symmetric Knight-shift tensor. The spectrum at 2.0 K 
has some broadening as manifested by smearing of the 
fine structure which appears at higher temperatures. As 
the temperature is decreased further, the line becomes 
broader and exhibits again a characteristic shape with 
two shoulders on both sides of the peak (except the small 
one near the zero shift which probably comes from a trace 
of nonmagnetic impurity phases). The broadening of the 
spectrum signals the onset of long-range magnetic order- 
ing of V 4+ spins giving rise to a finite internal magnetic 
field at the P site. The ordering is considered to be anti- 
ferromagnetic because the spectrum is broadened almost 
symmetrically about its center-of-gravity position in the 
paramagnetic state. As shown in the inset of Fig. 3, 
the FWHM of the spectrum exhibits a sudden increase 
around 2.0 K which gives a rough estimate of the Neel 
temperature Tn. Indeed, we determined Tn more pre- 
cisely from the T dependence of the nuclear spin-lattice 
relaxation rate 1/Ti at the P site to be 1.93(1) K where 
1/Ti is strongly peaked due to critical slowing down of 
the electronic spins (see Fig. 5). 

The NMR spectrum below Tn does not have a rect- 
angular shape which we usually observe for nonmag- 
netic nuclei in collinear antiferromagnets. The observed 
line shape, however, can be explained as the one in 



the two-sublattice antiferromagnet by taking account of 
1) the anisotropy of the hyperfine coupling resulting in 
the anisotropic Knight-shift tensor in the paramagnetic 
state, 2) incomplete cancellation of the transferred hy- 
perfine fields from neighboring V 4+ spins belonging to 
different sublattices due to the difference of the inter- 
atomic distances between V and P atoms, and 3) small 
tipping in the moment direction from the easy axis under 
the applied field of which effect is expressed by the par- 
allel and perpendicular susceptibilities. Details of the 
calculation and the simulation of the spectrum will be 
given in the appendix, and we only mention here that 
as shown in Fig. 11 the observed spectrum is well repro- 
duced using physically-reasonable parameters. 

From the values of the intrachain coupling J and the 
Neel temperature Tn, we can estimate the interchain cou- 
pling. The intoxhain coupling J' may be evaluated using 
the expression! 2 !]) 



J 



2J 



(3.2) 



Substituting Tn = 1.93 K and J/k-Q = 3.50 K into eq. 
(3.2), we obtain z\J'\/k-Q — 2.13 K. If we disregard for 
simplicity the difference between, the second and third 
nearest neighbor V-V distances,! 2 .!! z — 4 and we get the 
"average" interchain coupling \J'\/k& = 0.53 K (the sign 
cannot be specified). The ratio \J'\/J — 0.15 between 
the intra- and interchain exchange couplings is much 
larger than those of the other well-known examples of 
a ID magnet having J'/ J of order 10~ 2 or less. 

A careful analysis of the line shape in the paramagnetic 
state enables us an independent determination of the 
principal values of the Knight-shift tensor O We plotted 
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Fig. 6. Nuclear spin-lattice relaxation rate at the P site in a- 
VO(POs)2 as a function of H~ at two different temperatures. 
The dotted lines are the fit of the data to the formula 1/Ti = 
P + QH- 1 / 2 . 



Fig. 5. Temperature dependence of the nuclear spin-lattice re- 
laxation rate at the P site in o-VO(P03)2 measured at various 
magnetic fields. The inset shows the relaxation rate near Tn 
measured at 1.71 T. 



in Fig. 4 the T dependence of the two independent prin- 
cipal components of the Knight-shift tensor, K\\ and K±, 
which correspond to the shift for the magnetic-field di- 
rection parallel and perpendicular to the symmetry axis 
at the P site, respectively. The Knight-shift tensor de- 
pends strongly on temperature and takes a maximum at 

4.2 K as the bulk \- As shown in the inset of Fig. 4, 
both if || and K± scale with the bulk x above and below 
the susceptibility maximum. From the linear slopes of 
the if || and K± versus x plots, we determined the prin- 
cipal components of the hyperfine tensor at the P site 
to be A l{ = 5.5(1) kOe//; B and A ± = 3.3(1) kOe/// B . 
These values yield the isotropic and uniaxial components 
of the hyperfine coupling, ,4 iso = (A\\ + 2Ax)/3 = 4.1(1) 
kOe/^B and A ax = (A\\ - A x )/3 = 1.1(1) kOe/^ B , re- 
spectively. The hyperfine coupling at the P site is dom- 
inated by an isotropic transferred hyperfine field from 
the neighboring V 4+ spins and has a small anisotropic 
component. A ax is much larger than is expected from 
the classical dipolar field of the surrounding V 4+ spins 
(~ 0.18 kOe//i B ) and may be attributed to polarization 
of anisotropic p orbitals on the P atom. 

3.3 Nuclear spin-lattice relaxation 

Figure 5 shows the T dependence of the nuclear spin- 
lattice relaxation rate 1/Ti at the P site under various 
magnetic fields. The relaxation rate is strongly field 
dependent and is larger at lower fields. The field de- 
pendence becomes weaker as the temperature decreases, 
and 1/Ti becomes almost field-independent at 4.2 K. As 



to the T dependence, 1/Ti measured at different fields 
have some common characteristics. At high tempera- 
tures above about 50 K, 1/Ti is only weakly T depen- 
dent and decreases gradually on increasing temperature. 
On the other hand, the relaxation rate decreases rapidly 
below about 50 K except the rate at the field H = 4.34 
T showing a modest upturn in that temperature region. 
Note that the decrease of 1/Ti at low temperatures is 
more pronounced at lower fields. Below about 2.5 K, 
1/Ti exhibits a critical increase toward the long-range 
magnetic-ordering temperature as shown in the inset of 
Fig. 5. The maximum of 1/Ti is observed at 1.93 K 
which we determined as the Neel temperature of this 
compound. 

The strong magnetic-field (H) dependence of 1/Ti is 
likely to come from spin diffusion characteristic of a low- 
dimensional Heisenberg spin system. In one dimension 
the spectral density of the spin-spin correlations diverges 
as lo^ 1 / 2 toward lj — > which can be probed as H~ 1 / 2 
dependence of 1/Ti. Figure 6 shows examples of the H 
dependence of 1/Ti at the P site plotted as a function of 
H^ 1 / 2 . It is clear that 1/Ti at high temperatures obeys 
an H^ 1 / 2 law indicating diffusive behavior of the spin- 
spin correlations. The field dependence can be fitted to 
the form 

1/Ti = P + QH- 1 ' 2 (3.3) 

where P and Q are fitting constants. While the sec- 
ond term represents the contribution of spin diffusion 
near q — 0, the first term includes all possible field- 
independent contributions to 1/Ti. We analyzed the 
measured 1/Ti based on eq. ( fOj ) down to 4.2 K to 
parametrize the T and H dependence of 1/Ti, although 
the phenomenological theory of spin diffusion is justified 
in the limit T ^ J/fc B . Physical meanings of the ob- 
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Fig. 7. Temperature dependence of the field-independent relax- 
ation rate (l/Ti)o (open symbols) and the diffusive contribution 
(l/TiJdiff at 1.71 T (solid symbols). The solid lines are guides 
to the eyes. 



tained parameters will be examined in the next section. 
Returning to Fig. 6, the diffusive contribution propor- 
tional to iJ -1 / 2 dominates the nuclear-spin relaxation 
at 260 K. The slope in H~ x / 2 decreases monotonically 
with decreasing temperature while the field-independent 
part of 1/Ti grows up. The H~ x / 2 contribution vanishes 
almost entirely at 4.2 K, so that 1/Ti is governed by 
field-independent relaxation processes. Takigawa et al. 
reported a weak H dependence of 1/Ti in S^CuOs with 
a decreasing slope in the 1/Ti vs H~ x / 2 plot on cool- 
ingJiiP similar to the behavior of 1/Ti at the P site in 
a-VO(P03)2 at low temperatures. 

The T dependence of the field-independent relaxation 
rate (l/Xi)o = P is shown in Fig. 7. (l/Xi)o is almost 
T independent with the value of about 1.0 x 10 3 s _1 
above 50 K. Below about 50 K, on the other hand, a 
substantial increase of (l/Ti)o was observed. This is 
contrasted with the behavior of 1/Ti in that temperature 
region where 1/Ti decreases rapidly on cooling. The low- 
T decrease of 1/Ti therefore results from the decrease 
of diffusive contribution which overrides the increasing 
contribution of (l/7i)o- It is likely that the increase of 
(l/Ti)o from the high-T asymptotic value is caused by an 
enhancement of short-range AF correlations within the 
chain, because geometry at the P site allows intrachain 
AF spin fluctuations to contribute to the 31 P nuclear- 
spin relaxation. 

In Fig. 7 we also show for comparison the T de- 
pendence of diffusive contribution (1/Ti)diff = — 
(1/Ti) at 1.71 T. It is clear that, while govering 1/Ti 
at higher temperatures, (l/Ti)^ decreases rapidly be- 
low about 50 K. The decrease of (1/T 1 ) ( jiff synchronizes 
with the increase of (l/Ti)o, which strongly suggests 
some common origin for them. As mentioned above, the 
intrachain short-range AF correlation is a likely source 
for such T-dependent behavior of the nuclear-spin relax- 
ation. 



Fig. 8. Temperature dependence of the spin-diffusion constant 
D s . The dotted line is a value of D s in the classical limit. 



From the slope Q of the 1/Ti versus H~ x l 2 plot, we 
can estimate the spin-diffusion constant D s which de- 
termines the decay of the spin-spin correlation functions 
for small q as (Sq(t)S^_ q (0)) oc exp(— D s q 2 t) exp(— iuj c t) 
where ui« = gix^Hjh is the electron Larmor frc- 
quency.E-j If the hyperfine coupling is predominantly 
isotropic which is the case for the P site in a- 
VO(P03)2, the transverse component of the spectral 
density S xx (iu) = E, die™* ({S% (t) , S%(0)}) dom- 
inates the ^-dependent part of 1 jT\ . The nuclear spin- 
lattice relaxation rate (1/Ti)<jiff due to spin diffusion is 
then given as 



1 



A 2 k B T X 



spin 



Tij dm ti z g 2 ^^m^r c (3 ' 4) 

in one dimension. Here A is the hyperfine coupling con- 
stant in units of energy and Xspin is the spin suscep- 
tibility per magnetic atom. We determined D s using 
the values A iso for A, Xspin = {x ~ Xo)/N A and g de- 
termined from the analysis of the susceptibility and the 
Knight shift. The result is shown in Fig. 8 where D s 
is plotted as a function of temperature. D s is nearly 
T independent with the value of (6 — 8) x 10 11 s" 1 
above 20 K. It agrees well with the classical limitQ^ 
(7/^)^/2^5(5+ l)/3 = 5.7 x 10 11 s" 1 , which suggests 
that the spin dynamics is governed by the classical spin 
diffusion. On the other hand, D s becomes unusually 
large below about 10 K compared with the classical 
value. Such an anomalously large D s may be a signa- 
ture that spin diffusion no longer describes the intrinsic 
dynamics of the 5=1/2 HAFC at low temperatures. 

§4. Discussion 

Although a-VO(PC>3)2 undergoes a long-range AF 
transition at Tn = 1.93 K, the spin dynamics still pos- 
sesses a one-dimensional (ID) diffusive character at high 
temperatures as manifested by the strong H^ 1 / 2 depen- 
dence of 1/Ti. On the one hand, a weaker field depen- 
dence of 1/Ti below about 50 K implies low-energy spin 



6 



Jun Kikdchi et al. 



excitations being governed by different physics at lower 
temperatures. The change of the if-dependent behav- 
ior of 1/Ti around 50 K is nothing but an evidence for 
crossover between the two distinct temperature regimes 
of interest we mentioned in the Introduction, the param- 
agnetic regime T 3> J/k B at temperatures well above the 
intrachain exchange and the short-range-ordered (SRO) 
regime T ~ J/k B with strong intrachain AF correlations. 

In the paramagnetic regime above about 50 K, physi- 
cal quantities characterizing the dynamics such as 1/Ti, 
(l/Ti)o and D s are nearly T independent. 1/Ti is gov- 
erned by ID diffusive spin excitation near q = giv- 
ing rise to the H~ x / 2 dependence of 1/Ti at the P site. 
On the other hand, the field-independent relaxation rate 
(l/Ti)o remains small. The spin-diffusion constant D s in 
this temperature regime is in good agreement with the 
value expected in the classical limit, which suggests that 
the observed diffusive behavior is an intrinsic property 
of the spin system but not an effect of the coupling with 
the other degrees of freedom such as phonons£P 

It has been argued that the field-independent relax- 
ation rate in the paramagnetic diffusion regime comes 
from the longitudinal component of the spectral den- 
sity S zz (lu), of which divergence toward w — > is cutoff 
by interactions that break the conservation of the in- 
trachain uniform magnetization.ucP Here we estimate 
such a contribution to (l/Ti)o- The cutoff effect is taken 
into account phenomenologically by multiplying the ex- 
ponential decay to the longitudinal diffusive spin correla- 
tion function as (S*(t)S ! L q (0)) oc exn(— D s q 2 t) exp(— u c t) 
where uj c is the cutoff frequency .E3li 3 f On the assump- 
tion that the nuclear Larmor frequency is much smaller 
than u C) the relaxation rate (l/Ti) cuto ff due to the lon- 
gitudinal component with cutoff is given as 



A' k B T 



Xs 



pin 



cutoff & 3 2 MB VAsWc 

which in the limit T — > oo can be written as 



A' 2 5(5*+ 1) 



utoff 



n 2 3^ 



(4.1a) 



(4.1b) 



Here A' is the relevant coupling constant which in the 
present case is of dipolar origin. An important cutoff 
mechanism in a-VO(P03)2 is the interchain coupling 
much larger than the dipolar coupling between electronic 
spins. Taking uj c = k B J'/h « 7 x 10 10 s _1 and D s m 
7 x 10 11 s _1 , and using the dipolar coupling constant 
(~ 0.18 kOe/^B) for A 1 , we obtain (l/Ti) cutoff = 17 s" 1 
from eq. (4.1b). This is nearly two orders-of-magnitudc 
smaller than the observed (l/Ti)o- Even if we take 
smaller uj c due to the electron dipolar coupling (~ 0.05 
K), (l/Ti) cu toff is at most 54 s _1 and cannot explain 
the observed (1/Ti) . The above estimate of (l/Ti) cuto ff 
leads us to conclude that the field-independent relax- 
ation rate (l/Ti)o in the paramagentic regime is deter- 
mined by the mechanism other than spin diffusion, al- 
though the origin is not clear at present- It is noted that 
the exchange-narrowing limit of 1/TiEJ calculated us- 
ing parameters for a-VO(P03)2 as 1/Tioo = 2.6 x 10 3 
s" 1 cannot explain as well the value of (1/Ti) at high 



temperatures. 

Next we discuss the dynamics at low temperatures 
below about 50 K. An important observation in this 
SRO regime is a decreasing contribution of spin diffu- 
sion to 1 /Ti . Concurrently with a growth of intrachain 
AF correlations, diffusive excitations become less impor- 
tant as a possible channel for nuclear-spin relaxation, 
and at temperatures T J/k B where the AF alignment 
of the spins along the chain is almost established, the 
diffusive contribution becomes essentially absent. This 
accompanies an increase of D s as a fitting variable from 
the T-independent asymptotic value in the paramagnetic 
regime. The apparent increase of TJ S , however, seems to 
make no qu antitative sense because D s is estimated us- 
ing eq. (3.4) which is justified in the limit T 3> J/ks- 
It may rather suggests the classical spin-diffusion theory 
becoming inapplicable in the SRO regime: D s cannot be 
defined as a physical- meaningful spin-diffusion constant. 
This may result from a qualitative change of the intrin- 
sic dynamics of the S — 1/2 HAFC from diffusive to 
nondiffusive, possibly propagating ones as approaching 
T^J/k B . 

It should be noted that the breakdown of the clas- 
sical spin-diffusion theory due to SRO does not affect 
the qualitative behavior of (l/Ti)o- Since 1/Ti tends to 
be if -independent as the temperature decreases, it is ex- 
pected that the result for (l/Ti)o does not depend on the 
precise form of 1/Ti as a function of H . We may there- 
fore conclude that the low-T increase of (l/Ti)o is not 
an artifact but an intrinsic property of the system. Al- 
though we cannot exclude a possibility that the uniden- 
tified 77-independent contribution in the paramagnetic 
regime dominates (l/Ti)o in the SRO regime as well, it 
seems reasonable to consider that the increase of (l/Ti)o 
comes from a growing contribution of the AF excitation 
mode because 1 /Ti exhibits a divergent behavior toward 
T N . 

In summary, the low-frequency spin dynamics in a- 
VO(POs)2 is characterized by the two distinct tempera- 
ture regimes and the crossover behavior between them. 
In the high-T regime (T ^> J/k B ) 1/Ti is dominated 
by the ii-dependent diffusive contribution which is de- 
scribed quantitatively by the ID classical spin-diffusion 
theory. As the temperature approaches the intrachain 
coupling strength J/k B , the system goes into a differ- 
ent regime with an almost ii-independent 1/Ti. If we 
continue to analyze 1/Ti in terms of the classical spin- 
diffusion theory, we observe an apparent increase of the 
spin-diffusion constant from the classical limit. It rep- 
resents parametrically the decrease of diffusive contri- 
bution to 1/Ti, and may be interpreted as a sign of 
crossover to the low-T regime with nondiffusive, prop- 
agating spin excitations. 

Appendix: Analysis of the NMR Line Shape in 
the Antiferromagnetic State 

The NMR line shape observed in the antiferromag- 
netic (AF) state is unusual and is to be examined fur- 
ther. We present in this appendix the derivation of the 
resonance condition appropriate for this specific exam- 
ple and may be applicable if the hyperfine coupling at 
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the nuclear site under consideration has uniaxial sym- 
metry. We also give the calculation of the line shape 
using the derived resonance condition. From comparison 
of the calculated spectrum with the observed one, it is 
shown that the observed line shape below Tn is in consis- 
tency with collinear antiferromagnetic structure having 
two sublattices. The point of derivation is to take ac- 
count of an anisotropy of the hyperfme coupling, incom- 
plete cancellation of the transferred hyperfme fields from 
neighboring magnetic ions, and field-induced canting of 
sublattice magnetizations. 

The resonance frequency uj for nuclei with the nuclear 
gyromagnetic ratio 7 is generally given as 



w/~r=\Ho + H n \. 



(AT) 



Here if and if n are external and internal magnetic 
fields, respectively. Neglecting the classical dipolar field 
for the time being, if n is due to the hyperfme field and 
is given as 



£ 2, 



(A-2) 



where Sj is the electronic spin at the j-th site and Aij is 
the hyperfme tensor between the sites i and j. Sj fluc- 
tuates in time and hence if n in eq. (AT) should be the 
time average because we are dealing with the resonance 
frequency. For simplicity of notation we omit the brack- 
ets ( ) to show explicitly the time-averaging procedure 
and use Sj and if n as the time-averaged quantities in 
the following. 

In a-VO(PC>3)2, there are two neighboring V atoms in 
the distances 3.292 A and 3.392 A around the P site and 
the others are relatively far apart (> 4.613 A). Because a 
transferred hyperfme coupling is generally short-ranged, 
we may consider a contribution of these two V atoms to 
if n at the P site. Labelling these V sites as V(l) and 
V(2), and similarly the electronic spins as Si and S2, we 
can rewrite eq. (A-2) as 



H n = A 1 -S 1 + A 2 - S 2 



(A-3) 



for the P site in a-VO(PC>3)2. Here A\ and A 2 are the 
hyperfme tensors from V(l) and V(2) sites, respectively. 
Since V(l) and V(2) belong to the different AF chains, 
S\ and S2 can in principle be either parallel or antipar- 
allel in zero external field. 

If Si and S2 arc parallel, i.e., V(l) and V(2) be- 
long to the same sublattice, Si = S 2 and hence if n = 
(Ai + A 2 ) ■ Si. Because the sum A\ + A 2 corresponds 
to the hyperfme tensor in the paramagnetic state, there 
should be a transferred field of significant magnitude 
(a few kOe) at the P site. We expect in that case a 
usual rectangular shape of the spectrum which contra- 
dicts with the observation. We therefore assume antipar- 
allel alignment of Si and S 2 under zero external field 
which means that V(l) and V(2) belong to different sub- 
lattices. 

Under nonzero external field, the sublattice moment 
changes its direction and magnitude slightly from the 
zero-field values. If the external field is smaller than the 
spin-flop field, the effect may be expressed by the parallel 




easy axis 



Fig. 9. Schematic view of the canting due to the external field 
H of the electronic spins S 1 and S2 ■ The dotted line indicates 
the direction of the easy axis. Decomposition of Hq to the com- 
ponents parallel and perpendicular to the easy axis, Holl an d 
Hqj_, respectively, is also shown. 



and perpendicular susceptibilities, x\\ an d X-L, in the AF 
state. It is then convenient to decompose Si and S 2 into 
components parallel and perpendicular to the easy axis 
for the sublattice moment as Sj = Sj\\ + Sj± (j : 
and similarly the external field as Hq — H u +ifo_ 



= 1,2) 
(Fig. 



9). The following relations hold between the susceptibil- 
ities and the parallel and perpendicular components of 
Si, S 2 and if ; 



S i\\ 



2 X\\ H o\\ 



Si± + S 2 ± — 2x±H ^ 



(A-4a) 
(A-4b) 



Here Xll an d X-L are defined as the susceptibilities per 
mag neti c atom. Defining A = Ai+A 2 and 2 = A\ — A 2 , 
eq. (|A-3|) can be written as 



if n = X ±A ■ H 0± +x\\A- -Hoy 

+ l -[a ■ (Six - S 2± ) + a ■ (Si,, - 3 n )] (A-5) 

«Xx2.-ia" _L+X||2-i? || + 5-Si||. (A-6) 

On going from (A-5) to (A-6) we used the approxima- 
tions Siy w — S211 and Sij_ « S2j_. The first and second 
terms in eq. (A-6) correspond to the hyperfme field re- 
sulting from canting of Si and S2 from the easy axis 
under nonzero external field. The third term represents 
the hyperfme field due to the difference between A\ and 
A 2 . Putting a • Si || = H n0 in eq. (A-6) gives a more 
convenient form for H„: 



H n = X ±A ■ H ± + Xll -4 • if on + if no 



(A-7) 



As has been pointed out, if n o comes from uncancellation 
of the hyperfme fields from the two neighboring V atoms. 

In deriving the expression of the resonance condition 
from which the line shape is calculated, we have to take 
some orthogonal frame and specify the directions of ffo 
and if n o. It is reasonable to take the principal frame 
of the hyperfme tensor A at the P site in which A is 
diagonalized and the principal values are known from 
the experiment. Taking the z axis as the unique axis 
of the tensor A having axial symmetry, we denote the 
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Fig. 10. Principal frame used in the calculation of the resonance 
frequency. The dotted line with an arrow indicates the direction 
of the easy axis. The easy axis is taken to be in the zx plane. 



directions of the easy axis and the external field Hq by 
the polar and azimuth angles as shown in Fig. 10. The 
easy axis can be taken to lie in the zx plane so that 
(p a = without loss of generality because of the axial 
symmetry of A. For the field H n o we assume that it is 
in the direction of the easy axis for simplicity. This is 
equivalent to neglect the anisotropy of a. It is noted that 
with this assumption, the classical dipolar field may be 
included in H n Q so that H n o is regarded as the sum of 
classical dipolar and hyperfine fields, because the dipolar 
field is a likely source for the anisotropy field in the S — 
1/2 system and may be taken to be parallel to the easy 
axis. 

Our aim is to express the resonance frequency u as a 
function of the angles 0, tp and 9 n , and the field strength 
Ho = | -Ho | and H n0 = |H n o|. The angles and <p dis- 
tribute randomly in the polycrystalline sample, so that 
a specific line shape called powder pattern is observed. 
In the principal frame of A, Hq and H n are written 
explicitly as follows; 



A X 



A 



A = 



A», 



A 



An 



Hn 



H 



nO 




Ho sin cos ip N 
Ho sin sin ip 
Ho cos 



H n o sin 6 


-ffnO cos n 



(A- 



(A-8b) 



(A-8c) 



A|| and A± are the principal values of A which can be de- 
termined from the K-\ plots in the paramagnetic state. 
Hoy is just the projection of Ho onto H n o, and hence 
-ffoll an d Hq± can also be written down by components. 



Substituting eq. (A-7) into eq. (AT) and writing explic- 
itly the components in the principal frame of A using 
eqs. ( [A8| ), we get after lengthy calculations 

u? = uj x 2 + uj y 2 + uj 2 (A-9a) 



with 

Ux/l = [{1 + A±(c x 2 x\\ + c z 2 x±)}b x 

+A±(x\\ - X±)cxC z b z ]H Q + c x H a0 , 
uyh = (1 + A±x±)b y Ho, (A-9b) 
w z /7 = [{1 + A ± (c z 2 X \\ + c x 2 x ±)}b z 

+A||(x|| - x±)cxC z b x ]H + c z H n0 . 

Note that x\\ (x±) i s not the susceptibility for the z (x, y) 
direction(s) in the principal frame of A but the one paral- 
lel (perpendicular) to the easy axis. Neglecting the terms 
of 0(A 2 x 2 ), and using the relations b x 2 + b y 2 + b 2 = 1 
and c x 2 + c 2 = 1, eq. (|A-9|) can be rewritten as 



Lu 2 /~f 2 = H Q 2 + H n0 2 + 2{b x c x + b z c z )H H n0 

+ 2[Aj_(c x 2 x\\ + c 2 x±)b x 2 + Aj_x±b y 2 

+A ll (c z 2 x\\+c x 2 X ±)b z 2 

+(A|| +A±)(x\\ - X±)b x b z c x c z ]H 2 

+ 2[A ± (c x 2 x\\ + c z 2 x ±)b x c x (A- 10) 

+A||(c z 2 X|| + c x 2 x±)b z c z 

+A||(x|| - x±)cz 2 b x c x 

+A±(X\\ - X±)cx 2 b z c z ]H H n o. 

Equation (ATO) is valid for arbitrary strength of H n o as 
far as the conditions Ax <§C 1 are satisfied. 

If H n o is much larger than Ax, eq. (ATO) may be 
simplified as 

uj 2 h 2 = Ho 2 + H n0 2 + 2{b x c x + b z c z )H H n0 (A- 11a) 

which, by denoting the angle between Hq and H n o as 
(3, reduces to 



u 2 h 2 = Ho 2 + H n0 2 + 2H H n0 cos/3. 



(ATlb) 



Random distribution of [3 in the polycrystalline sample 
yields a usual rectangular shape of the spectrum.ES 

If H n o is much smaller than Ho but is comparable with 
Ax which is the case here, we may neglect in eq. (ATO) 
the terms proportional to H n o 2 or AxH n o- Hence we 
arrive at the final expression of the resonance frequency; 

w/7 = H + [A ± (c x 2 x\\ + c 2 xi_)b x 2 + A ± x±b y 2 

+A ll (c z 2 X \\+c x 2 x±)b z 2 (AT2) 

+(A|| + Ax)(X|| - X±)bxb z c x c z ]H a 

+ (b x c x + b z c z )H n0 , 

or writing explicitly the dependence on the angles as 

w/7 = H + [Aj_(x\\ sin 2 n + x± cos 2 n ) sin 2 ^cos 2 (p 



+ Aj_xj_ sin 2 sin 2 ip 

+Aii (x|i cos 2 D + x± sin 2 n ) cos 2 



(ATS) 



+ -(A|, + A ± )(x|| -X±)sin20 n sm20 cos (p]H 
+ (sin n sin cos <p + cos n cos 0)H nO . 
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Fig. 11. Comparison of the calculated powder patterns with the 
observed spectrum at 1.65 K. The solid circles represent the ex- 
perimental data. The dashed line is the calculated spectrum 
f(Ho) with no inhomogeneous broadening. The parameters 
used in the calculation are uj = 48.5 MHz, An = 5.5 kOe//^B, 
A ± = 3.3 kOe/^B, H n0 = 310 Oe, 9 n = 33°, x\\ = 0.026 
cmu/mol and \± = 0.031 emu/mol. The spectrum shown by 
the solid line is calculated by taking account of inhomogeneous 
distribution of the internal field H n Q. The effect is introduced 
by convoluting lorentzian with the FWHM of 55 Oe. 
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Sj fluctuates in time and hence H n in eq. (AT) should be 



the time average because we are dealing with the resonance 
frequency. For simplicity of notation we omit the brackets ( ) 
to show explicitly the time-averaging procedure and use Sj 
and H n as the time-averaged quantities. 



For the field-sweep measurement the resonance fre- 
quency ui is fixed while the external field Hq is varied. 
It is therefore necessary to solve eq. (A- 13) with respect 
to Hq , the result of which is not given here explicitly be- 
cause it is rather straightforward. Then we calculate the 
spectrum /(Hq) oc AN/ AHq by counting the number of 
nuclei AN having the resonance field between Ho and 
Hq + AHq. More than million different (6, if) points rep- 
resenting the random distribution of Hq were taken to 
calculate the resonance field. The unknowns H n o, 6 nj 
and x± m (A- 13) are treated as parameters to reproduce 
the observed spectrum. As shown in Fig. 11, an ap- 
propriate choice of the parameter values reproduces the 
observed two-shoulder structure of the spectrum in the 
AF state. We may therefore conclude that the magnetic 
ordering in a-VO(PC>3)2 is not unusual but is rather con- 
ventional. 
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